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Abstract 

We analyse MINBU distribution of 2 dimensional quantum gravity. New 
data of R 2 -gravity by the Monte Carlo simulation and its theoretical analysis 
by the semiclassical approach are presented. The cross-over phenomenon takes 
place at some size of the baby universe where the randomness competes with 
the smoothing force of i? 2 -term. The dependence on the central charge c m 
and on the i? 2 -coupling are explained for the ordinary 2d quantum gravity 
and for i? 2 -gravity. The i? 2 -Liouville solution plays the central role in the 
semiclassical analysis. A total derivative term (surface term) and the infrared 
regularization play important roles . The surface topology is that of a sphere 



1 Introduction 

The quantum effects of the 2 dimensional (2d) gravitational theories are recently 
measured numerically in the computer simulation with high statistics. In particular 
the data for the entropy exponent (string susceptibility) in 2d quantum gravity(QG) 
is the same as the known exact result within a relative precision of 0(1CU 3 ). It is due 
to the developement of the simulation technique in the dynamical triangulationJT], 
||] and the findings of new observables in QG such as MINBU distribution 0, ||, ||. 

The data analysis is done by a rather orthodox approach, i.e., the semiclassical 
approximation. It has recently been applied to 2d _R 2 -gravity and the simulation 
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data of < / d 2 x^fg~R 2 > and its cross-over phenomenon are successfully explained |7|]. 
We list the merits of this approach. 



1. The semiclassical treatment is, at present, the unique field-theoretical ap- 
proach which can analyse the mysterious region (25 >)c m > 1. The conformal 
field theory gives a meaningful result only for some limitted regions of c m . The 
Matrix model is in the similar situation. 

2. Comparison with the ordinary quantization is transparent because the ordi- 
nary renormalizable field theories ,such as QED and QCD, are quantized es- 
sentially in the semiclassical way. In particular, the renormalization properties 
of (2d) QG are expected to be clarified in the semiclassical approch||. 

3. This approach can be used for the higher- dimensional QG such as 3d and 4d 
QG. 

The approach is perturbative, therefore choosing the most appropriate vacuum under 
the global constraints (such as the area constraint and the topology constraint) is 
crucial in the proper evaluation. We explain it in Sect. 3. 

We add i? 2 -term to the ordinary 2d gravity for the following reasons. ( We call 
the ordinary 2d gravity Liouville gravity in contrast with R 2 -gravity for the added 
one. ) 

1. For the positive coupling, the term plays the role of suppressing the high 
curvature and making the surface smooth. For the negative one, the high 
curvature is energetically favoured and making the surface rough. Therefore 
we can expect a richer phase structure of the surface configuration. 

2. The term is higher- derivative (<9 4 ), therefore it regularizes the ultra-violet be- 
haviour so good|J. In fact the theory is renormalizable ||. 

3. The Einstein term (R-term) is topological in 2 dimension. It does not have 
a local mode. The simplest interaction which is purely geometrical and has 
local modes is i? 2 -term. 

4. In the lattice gravity, i? 2 -term is considered as one of natural irrelevent terms 
in the continuous limit ||10||. 



The < / d 2 x^fg~R 2 > simulation data for i? 2 -gravity was presented by |TT] and the 
cross-over phenomenon was clearly found. We present here MINBU distribution 
data. 



2 Lattice Simulation of 2D Quantum R 2 -Gravity 
and MINBU Distribution 

The distribution of baby universe (BU) is one of important observables in the lattice 
gravity @, |], ||. It was originally introduced to measure the entropy exponent (string 
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susceptibility) efficiently. Fig.l shows the configuration of a BU with an area B 
(variable) from the mother universe with an area A (fixed). 



Fig.l MINBU configuration 

The 'neck' of Fig.l is composed of three links which is the minimum loop in the dy- 
namically triangulated surface. The configuration is called the minimum neck baby 
universe (MINBU). MINBU distribution for the Liouville-gravity and its matter- 
coupled case were already measured |5], ||, [13 



First we explain briefly our lattice model of i? 2 -gravity. The surface is regularized 
by the triangulation. The number of vertices ,where some links (edges of triangles) 
meet, is Nq. The number of links at the i-th vertex (i = 1,2 • • • , No) is The 
number of triangles (N 2 ) is related to No as iV2 = 2N — 4 for the sphere topology. 
The discretized model is then described by 

Sl = -Pl*£E& ^ = -487T 2 ^E4 + const ' W 

where (3l is the R 2 -coupling constant of the lattice model. We do measurement 
for p L = 0, 50, 100, 200, 300, -20, -50 . We present the MINBU dstribution of R 2 - 
gravity with no matter field (pure R 2 -gravity) in Fig. 2 and 3 for (3l > and for 
(3l < respectively. The total number of triangles is N% = 5000. For the detail see 



Fig.2 MINBU distribution for (3 L > 0, Pure i^-gravity. 
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Fig. 3 MINBU distribution for j3 L < 0. Pure i^-gravity. 

As for positive (3l (Fig. 2), we see clearly the transition point P , for each curve, 
at which the distribution qualitatively changes. For the region P = B/A > Pq, the 
birth probability decreases as the size of BU increases. For the region P < Pq, the 
birth probability increases as the size of BU increases. The value of the transition 
point P depends on (3 and increases as (3 increases. As for negative j3 (Fig. 3), the 
slope of the curve tends to be sharp as increases at least for the region P < Pi. 
The transition point Pi is not so clear as Fig. 2. 

In Sect 4.2 we interpret these data theoretically using the semiclassical approach 
explained in Sect 3. 



3 Semiclassical Approach 

We analyse the simulation data by the semiclassical approach. The R 2 - gravity 
interacting with c m -components scalar matter fields is described by 

S = Jd 2 x^(±R-(3R 2 -fi-lEi=ida®i-g ab -d b ® l ) , ( 0,6=1,2) ,(2) 

where G is the gravitaional coupling constant, \x is the cosmological constant , (3 is 
the coupling strength for i? 2 -term and $ is the c m - components scalar matter fields. 
The signature is Euclidean. The partition function , under the fixed area condition 



A — J d Xy/g and with the conformal-flat gauge g a b = e v S a b , is written as |L4 



Z[A] = J m^{ eX p±S} 5(J d 2 x^g - A) = e*p±(^^ - fiA) x Z\A\ , 

Z[A] = JVip e+i 50 ^ 6{f d 2 x e* - A) , (3) 
S^} = J d 2 x (^d 2 ^ - (3 e-^(d 2 ^) 2 + ^d a ( V d aV ) ) , I = i i-(26-c m ) (4) 

where h is the number of handles Q. Vgc is the gauge volume due to the general 
coordinate invariance. £ is a free parameter. The total derivative term generally 
appears when integrating out the anomaly equation 5Si n d[(p]/5(p = ^d 2 <p . This term 



The sign for the action is different from the usual convention as seen in (fj). 
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turns out to be very important. Q We consider the manifold of a fixed topology of 
the sphere ,h = and the case 7 > (c m < 26). Ti is Planck constant. 

Z[A] is rewritten as, after the Laplace transformation and the inverse Laplace 
one, 

Z[A] =SffV<pexpi[ S [<p] - X(fd 2 xe* - A)} 
S x y} = S [if]-Xfd 2 x 



Jd 2 x (^dV-/?e-*(dV) 2 + 4;d a (<pd a <p) -Xe 



27 



(5) 



where the A-integral should be carried out along an appropriate contour parallel to 
the imaginary axis in the complex A-plane. Note that the 5-function constraint in 
(|3]) is substituted by the A-integral. The leading order configuration is given by the 
stationary minimum. 



Sip 



-6 2 <p + P{e-*(d 2 <p) 2 - 2d 2 (e-*d 2 ip)} - Xe" 
7 







d 

dX 



(XA + S x [ Vc \) 



Z[A] » l -ex V \{X C A + S Xc [ Vc }} = l -ex V -if 



(6) 



Generally this approximation is valid for a large system. In the present case, the 
system size is proportional to — = 26 ~ 2 Cm . We expect the approximation is valid 
except the region: c m ~ 26. 

The solution (p c and A c , which describes the positive-constant curvature solution 
and is continuous at (3 = 0, are given by [J/J 



Mr) = {^(1 + -j?} 



(x 1 ) 2 + (x 2 ) 2 



An 1 

a c = — {w + 1 - \ w 2 + 1 - 2£w } 
w v 



w = 167r/3'7 , 

X A W , 

lX c A = — (a c 

167T 



A 



(7) 



where £ must satisfy — 1 < £ < +1 for the realness of a c . (x 1 ^ 2 ) are the flat 
(plane) coordinates. The partition function at the classical level is given by 



Y e Jf = In Z[A]\ h o = X C A+ (1 



C(A) 



\h° 

87r(2+g) 

7 



-^ W + C(A) 



+ ^{ ln(L 2 /A) - 1 } + 0(A/L 2 



L 2 



> 1 



.4 

where L is the infrared cut-off (r 2 < L 2 ) introduced for the divergent volume integral 
of the total derivative term. Note that C(A) does not depend on (3 (or w) but 



2 The uniqueness of this term, among all possible total derivatives, is shown in Q. 

3 In this section only, we explicitly write Ti (Planck constant) in order to show the perturbation 
structure clearly. 
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on c m (or 7) and A. Furthermore C(A) has an arbitrary constant of the form 
(87r£/7) x (const) due to the freedom of the choice of the regularization parameter: 
L — > (const) x L. This arbitrary constant turns out to be important. 

For the case (3 = , the theory is ordinary 2d gravity and we call it Liouville 
gravity in contrast with i? 2 -gravity for (3^0. For the case c m = , the theory is 
called the pure gravity in contrast with the matter-coupled gravity c m 7^ 0. 



4 Semiclassical Analysis of MINBU Distribution 

First we explain the free parameter £. Recent analysis of the present theory at 
the (1-loop) quantum level has revealed that it is conformal (the renormalization 
group beta functions=0) for w > 1 when we take £ = 1 ||. Therefore the value 
£ = 1 has some meaning purely within the theory. The validity of this choice is also 
confirmed from a different approach, that is, the comparison of the special case /3(or 
w)= (Liouville gravity) of the present result with the corresponding result from 
the conformal field theory (KPZ result) ||. The asymptotic behaviour of Z[A]| a o at 
w = is given, from (M), as 
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= ezp{^(3 -{) + (1 + O-'n 1 -^ + M,4j 

w=o 7 7 2 7 A 

8££ 26- cm f 

A t x const = A s s x const 



as A -> +00 . (9) 

On the other hand, the KPZ result is 

Z KPZ L4]~^~ 3 , % = -^{cm - 25 - ^(25 - c m )(l - c m )} + 2 . (10) 

In order for our result to coincide with the KPZ result in the 'classical limit' c m — > 
-00 : Z KPZ [A] ~ A + l Cm , we must take 

e = i , (ii) 

in (0). In the following of this text we take this value. 

The asymptotic behaviour of the present semiclassical result for the Liouville 
gravity is, taking £ = 1 in (0), 

Z[A] ~ A~^ xA- 1 , A — +00 , (12) 

where the additional factor A -1 comes from the A-integral in the expression of 
Z[A], ([5])[|1^]. Now we compare the KPZ result and the semiclassical result in the 
normalized form. 

Z KPZr A l _ f l^i ^ 4 7 s(c m )-7 s (cm=0) 

norml j - r^KPZ^^ 

ls{c m ) - 7 s (c m = 0) = J^Wm + 5 - -y/(25 - c m )(l - c m )} , (13) 

Z^[A] = -$£-~A+r . 

Z L4llcm=0 



In the numerical evaluation, we take £ = 0.99 for the practical reason. 
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We can numerically confirm that the semiclassical result, and the KPZ result, 
7s (c m ) — 7 S (c m = 0) , have very similar behaviour for the region c m < 1[|1^ . 



Now we go back to the general value of (3. The birth-probability of the baby 
universe with area B(0 < B < A/2) from the mother universe with the total area 
A is given by pi] 



, m 3 {A - B + a 2 ) (B + a 2 )Z[5 + a 2 ]Z[A - B + a 2 } 

Ha{B) = A 2 x Z[A] 

3(1 - p)pZ\pA]Z[{l - p)A] 

In (^^1) « In (1 - p)p + In Z[pA] + In Z[(l - p)A] - In . 
3 

B 1 
p=— , < p < - . 
A y 2 

We apply the result of Z[A] in Sect. 3 to the above expressions. 



(14) 



4.1 c TO - dependence 

First we present the semiclassical prediction for Liouville gravity(/5 = 0). The result 
(H) for the case (3 = gives ,taking £ = 1, 

1 L 2 

7 Zn Z[rA] = 87r(ln n + 1) + 8vr ln( -) . (15) 

T J\ 

Then the MINBU distribution normalized by the pure garvity (c m = 0) is obtained 



as 



- {p(l — p)} e x exp {— x A} 



n A (B)\ Cm=0 ^ r L 12 

L 2 

A = -2(ln 7T + 1) -2 In— , (16) 

where A can be regarded as the free real parameter due to the arbitrariness of the 
infrared regularization parameter L. We know from the result (|16"D that the MINBU 
distribution lines for different c m 's cross at the single point p = p* given by 

p*(l —p*) = exp{-i A} , p*< 1 - . (17) 
Fig. 4 shows three typical cases of p* . 
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Fig.4 Three typical cases of the solution of ([17]). 

The choice of A is important to fit the theoretical curve (|T6| ) with the data. We show 
the behaviour of ( |l6l) for the three cases: 1) exp(— -A) <C | ,Near Point O, Fig. 5a 
; 2) exp(-±A) > \ ,Above Point A ,Fig.5b ; 3) exp(-±A) = § - ,Near Point A 
, Fig. 5c. 



Fig. 5a MINBU distribution for Liouville gravity, A = 8 



Fig. 5b MINBU distribution for Liouville gravity, A = 1 
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Fig. 5c MINBU distribution for Liouville gravity, A = 3 



Fig. 5a well fits with the known result of the computer simulation 0, |T3|. This 
result shows the importance of the infrared regularization. 



4.2 /^-dependence 

We consider the pure gravity(c m = 0). We plot MINBU dstribution, In Ua{B), as 
the function of p ( 0.001 < p < 0.1 ) for various cases of /3' = /3/A (£ = 0.99). Fig. 6a 
and 6b show that for /?' > and f3' < respectively. 



Fig.6a MINBU distribution for ft > 0. £ = 0.99, c m = 0. 



Fig.6b MINBU distribution for < 0. f = 0.99, c m = 0. 
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The above results of Fig.6a and Fig. 6b qualitatively coincide with those of Fig.2 
and Fig.3, respectively. 



We list the asymptotic behaviour of In ua{B) for the general £ and c m in Table 

1. 



Phase 


(C) <p< 1) 


(B) \w\ <p 


(A) < p < iw(~ 1) 


a~(pA) 


Mi + ¥S 


4 7 r(l + £){l- 1 ^ 


47r(i+e)p x 


In n A (B) 


(l-^)lnp 

47T UI 

7 P 

+SmallTerm 


(l-^)lnp 
+SmallTerm 


{1-Mi=ii }lnp 

4 ^ 1+ « In «, 

+0(Z) 
+ Small Term 



Table 1 Asymp. behaviour of MINBU distribution, (|14]) , 
for general c m and £. > 0, w = 16ir(3'j, 7 = OR 487r > 0, p = 4 , 
0<p<l, |w| ~ 1, SmallTerm = const + 0(wp) + 0(p). 



We characterize each phase in Table 1 as follows. 

(A) 0<p<w: Smoothly Creased Surface (] 

The smoothing term, R 2 , dominates the main configuration and the surface is 
smooth. The left part P < Pq(w) for each curve (w) in Fig. 6a corresponds to this 
phase. The small BU is harder to be born because it needs high-curvature locally. 
The large BU is energetically preferable to be born. The area constraint is not 
effective in this phase. The characteristic scale is (3. 

(B) |iy| <C p: Fractal Surface 

The randomness dominates the configuration. The size of BU is so enough large 
that the i? 2 -term is not effective. The area constraint is neither effective. There is 
no characteristic scale. The right part P > Po(w) for each curve (w) in Fig.6a and 
the right part P > Pi(w) for each curve (w) in Fig. 6b correspond to this phase. 
The MINBU distribution is mainly determined by the random distribution of the 



surface configuration |Rj| 



(C) < p — w: Rough Surface | 

Due to the large negative value of _R 2 -coupling, the configuration with the large 
curvature is energetically preferable on the one hand, it is strongly influenced by 
the area constraint on the other hand. Therefore the large BU is much harder to 
be born than (B) because it has a small curvature and a large area. The small BU 



5 In we called it Free Creased Surface because this is the phase where the free kinetic term 
(i? 2 -term) dominates. 

6 In we called it Strongly Tensed Perfect Sphere because the surface tension is negatively large 
and the shape of the whole surface is near a sphere. At the same time the surface tend to become 
sharp-pointed because it increases the curvature. We call the surface under this circumstace, simply, 
Rough Surface. 
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is much easier to be born than (B) because it has a large curvature and a small 
area. The left part P < P\{w) for each curve (w) in Fig.6b corresponds to this 
phase. The characteristic scale is the total area A. 



We see the phase structure of Table 1 is the same as that of by the 
substitution of w by w/p . Although both simulations measure the same surface 
property, the cross-over phenomenon, however, appears differently. In ]7) the 
physical quantity < / d 2 x^/gR 2 > is taken to see the surface property. The 
cross-over can be seen only by measuring for a range of w and the transition point 
is given by a certain value |iu*| m 1. This is contrasting with the present case. The 
cross-over can be seen for any w. The transition is seen at the point p* ,in the 
MINBU distribution, given by \w\/p* ~ 1. We understand as follows. The MINBU 
distribution measures the surface at many different 'scales' B, whereas the 
quantity < / d 2 Xy/gR 2 > measures the surface at a fixed ' scale '( B\ (or pi) in the 
MINBU terminology). 



4.3 General Case 

We consider the general case of c m and f3. This general case is not yet measured by 
the Monte Carlo simulation. We present the semiclassical prediction. The analysis 
so far shows the normalization (([13) and fll6|) ) and the choice of an arbitrary 
constant due to the infrared regularization ([16]) are important for the quantitative 
adjustment. Here, however, we are content with the qualitative behaviour. We 
donot do the normalization and we ignore the In ^- term in the evaluation of this 
subsection. 



(1) c m -dependence 

We stereographically show MINBU distributions for the range: 

0.001 < p < 0.2, -24 < c m < +24 , in Fig.7a(/3' = 0) , Fig.7b(/3' = +10" 4 ) and 

Fig.7c(/3' = -10" 5 ). 



Fig. 7a MINBU dstribution for 0.001 < p < 0.2, -30 < c m < +24 . 

/3> = 0^ = 0.99. 
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Fig.7b MINBU dstribution for 0.001 < p < 0.2, -30 < c m < +24 . 

P> = +1Q- 4 ,^ = 0.99. 



Fig. 7c MINBU dstribution for 0.001 < p < 0.2, -30 < c m < +24 . 

P> = -10~ 5 ,£ = 0.99. 

No 'ridge' appears in Fig. 7a. From this, we see matter fields affect the surface 
dynamics homogeneously at all scales. (This result is natural because the matter 
coupling constand c m does not have the scale dimension.) The slope along the 
p-axis continuously decreases as c m increases. In Fig. 7b, a ridge runs from a low p 
to a high p as c m increases. In Fig. 7c, a 'hollow' runs from a high p to a low p as c„ 
increases. The ridge and the hollow correspond to the series of the cross-over 
points. In both Fig. 7b and Fig. 7c, the cross-over becomes dimmer as c m increases 
and becomes sharper as c m decreases. 

(2) /5-dependence 

We stereographically show MINBU distributions for the range: 

0.001 < p < 0.2, -10~ 5 < P' < +10~ 4 , in Fig.8a(c m = 0) , Fig.8b(c m = +10) and 
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Fig.8c(c m = -10). 



Fig.8a MINBU dstribution for 0.001 < p < 0.2, -10~ 5 < (3' < +10" 4 . 

c m = 0,e = 0.99. 



Fig.8b MINBU dstribution for 0.001 < p < 0.2, -10~ 5 < /3' < + KT 4 . 

c m = +10,e = 0.99. 



Fig.8c MINBU dstribution for 0.001 < p < 0.2, -10" 5 < 0' < +10~ 4 . 

c m = -10,e = 0.99. 

The Fig. 8a corresponds to the stereographic display of Fig.6a and 6b. In each of 
Fig.8a-c, a ridge appears for f3' > . For (3' < , a tower appears instead of a 
ridge. For a large positive c m ( matter dominated region, c m = 10 in Fig.8b) the 
undulation of the MINBU dstribution surface [] is small(the cross-over is dim), 
whereas it is large (the cross-over is sharp) for a large negative c m (matter 
anti-dominated region, c m = — 10 in Fig. 8c). 

7 Do not confuse it with the 2d manifold which the present model of gravity represents. 
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5 Discussion and Conclusion 



In the (2d) QG,at present, there exists no simple way to find good physical 
observables. They have been found by 'try and error'. MINBU is one of good 
observables to measure the surface property. Quite recently a new observable ,the 
'electric resistivity' of the surface, is proposed by (T^|. By measuring the 



observable for the matter-coupled Liouville gravity, they observe a cross-over ,near 
c m = 1 ,from the surface where a complex-structure is well-defined to the surface 
where it is not well-defined. The analysis of the new obserbable, from the 
standpoint of the present approach, is important. 

There are some straightforward but important applications of the present 
analysis : 1) higher-genus case, 2) the case with other higher-derivative terms such 
as R 3 and Vi? • Vi? , 3) the quantum effect. As for 2) ,references [17] and |18| have 
already obtained the Monte Carlo data. 

We have presented the numerical result of MINBU and its theoretical 
explanation using the semiclassical approximation. The surface properties are 
characterized. It is confirmed that the present lowest approximation is very 
efficient to analyse 2d quantum gravity, at least, qualitatively. 

Finally we expect other new observables will be found and many Monte Carlo 
measurements will be done including 3 and 4 dimensional cases, next a few years. 
The interplay between the measurement by the computer simulation and the 
theoretical interpretation will become important more and more. We believe this 
process will lead to the right understanding of the (Euclidean) quantum gravity. 
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gure Captions 

• Fig.l MINBU configuration. 

• Fig. 2 MINBU distribution for (3 L > 0, Pure # 2 -gravity. 

• Fig. 3 MINBU distribution for (3 L < 0. Pure i^-gravity. 

• Fig.4 Three typical cases of the solution of ([17]). 

• Fig. 5a MINBU distribution for Liouville gravity, A = 8. 

• Fig.5b MINBU distribution for Liouville gravity, A = 1. 

• Fig. 5c MINBU distribution for Liouville gravity, A = 3. 

• Fig.6a MINBU distribution for p > 0. f = 0.99, c m = 0. 

• Fig.6b MINBU distribution for /3' < 0. £ = 0.99, c m = 0. 

• Fig. 7a MINBU dstribution for 0.001 < p < 0.2, -30 < c m < +24 . 
13' = o,f = 0.99. 

• Fig. 7b MINBU dstribution for 0.001 < p < 0.2, -30 < c m < +24 . 

(3' = +10" 4 ,^ = 0.99. 

• Fig. 7c MINBU dstribution for 0.001 < p < 0.2, -30 < c m < +24 . 
13' = -10" 5 ,£ = 0.99. 

• Fig.8a MINBU dstribution for 0.001 < p < 0.2, -10" 5 < f3' < +10" 
c m = 0,e = 0.99. 

• Fig.8b MINBU dstribution for 0.001 < p < 0.2, -10" 5 < (3' < +10" 
c m = +10,^ = 0.99. 

• Fig.8c MINBU dstribution for 0.001 < p < 0.2, -10" 5 < (3' < +10" 
c m = -10,^ = 0.99. 
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